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Final  Report  for  the  AFRL  grant  #  FA8650-05-1-6641 

Experiments,  Computation,  and  Modeling  for  Temperature  Dependence  of 
Absorption,  Scattering,  Reflection,  Transmission,  and  Index  of  Refraction  of 
Optical  Radiation  in  Biological  Tissues 

The  objective  of  this  research  was  to  systematically  measure  the  temperature 
dependence  of  optical  properties  of  various  bovine  eye  tissues  using  different 
experimental  and  computational  techniques  at  elevated  temperatures  (above  body 
temperature).  An  attempt  was  made  to  measure  the  temperature  dependence  of 
absorption  and  scattering  coefficients  of  bovine  ocular  tissue  in  the  temperature  range  of 
24  to  40°C  at  980  mn.  The  rate  of  change  of  the  refractive  index  with  respect  to 
temperature  was  perfonned  at  632.8  mn  for  bovine  aqueous  and  vitreous  humor  in  the 
temperature  range  of  37  to  41  °C. 


Temperature  Dependence  of  Absorption  and  Scattering  Coefficient 
Experimental  Details 

The  schematic  of  the  experimental  setup  for  measuring  the  temperature  dependent 
total  diffuse  reflectance  and  total  transmittance  of  bovine  ocular  tissue  is  shown  in  Fig.  1 . 
Two  identical  integrating  spheres  (Newport  70674)  each  of  diameter  150  mm  with  port 
apertures  of  1 1  mm  in  diameter  were  utilized.  A  laser  beam  was  split  by  a  30-70  beam 
splitter  (BS)  which  directed  30  percent  of  the  light  to  the  reference  sphere  and  the 
remaining  70  percent  of  the  light  toward  entrance  port  I  of  the  sample  sphere.  The  tissue 
under  study  was  placed  in  a  temperature-controlled  sample  holder  whose  details  are 
shown  in  Fig.  2.  The  sample  holder  was  made  from  6.44  mm  thick  aluminum  and  had 
two  1  mm  thick  glass  windows.  The  aperture  of  the  sample  chamber  was  approximately 
16  mm  in  diameter.  The  sample  holder  had  two  heaters  (Omega  SRFG101)  attached  to 
one  side.  The  temperature  of  the  sample  was  elevated  by  controlling  the  voltage  across 
the  heating  elements  with  a  variable  voltage  supply  (Powerstat  3PN1 16B)  and  measuring 
the  temperature  with  a  digital  thermometer  (Hart  Scientific  1006).  For  total  transmission 
measurements,  the  sample  holder  was  mounted  to  port  I  of  the  sample  sphere  while  a  port 
II  of  the  sample  sphere  was  covered  with  a  cap  with  a  reflective  surface  identical  to  that 
of  the  integrating  spheres.  For  total  diffused  reflectance  measurements,  the  sample 
holder  was  mounted  to  port  II  of  the  sample  sphere  while  a  port  I  was  open.  The  sample 
holder  was  thermally  insulated  from  the  sample  sphere  by  a  1  mm  thick  polycarbonate 
spacer. 


The  measurements  were  performed  on  bovine  vitreous  humor  at  980  mn  in  the 
temperature  range  of  21  -  40°C.  The  980  mn  wavelength  was  supplied  by  pig-tailed 
diode  lasers  (Thorlabs  model  L9805E2P5)  coupled  to  a  collimator  (Thorlabs  model 
F220FC-C)  and  controlled  by  a  Thorlabs  (model  LDC  2000)  laser  driver.  The  laser  beam 
diameter  at  l/e~  of  the  peak  intensity  was  1.97  mm  and  beam  divergence  was  about 
0.048°.  The  average  output  power  was  kept  at  about  0.25  mW  for  all  optical 
measurements. 
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The  reflected  and  transmitted  light  intensities  were  detected  by  two  identical 
infrared  detectors  (Judson  J16-D)  attached  to  the  two  measuring  ports  of  the  integrating 
spheres.  The  detector  signal  was  sent  to  identical  preamplifiers  (Judson  PA-9)  powered 
by  a  common  power  supply  (Pasco  8000).  Signals  from  the  detectors  were  measured  and 
recorded  by  a  digital  oscilloscope  (Tektronix  TDS3054B).  The  measured  light  intensities 
were  then  utilized  to  determine  the  total  diffuse  reflectance  Rd  and  total  transmittance  Tt 
by  the  following  expressions: 


and 


R, 


X,-Cxl 
2-33  (Z.-Cj 


z  = 


XT  CXT 

2.33  (Zr-CZT) 


(1) 

(2) 


where  Xr  is  the  measurement  at  the  sample  sphere  with  the  sample  at  port  II,  Zr  is  the 
measurement  at  the  reference  sphere,  Cxr  is  the  correction  factor  for  stray  light  measured 
at  the  sample  sphere  with  no  sample  at  port  I  and  port  II  uncovered  and  the  laser  on, 
while  Czr  is  the  correction  for  the  reference  sphere  with  the  laser  off.  For  the  total 
transmittance  measurements,  Xj  is  the  measurement  at  the  sample  sphere  with  the  sample 
at  port  I,  Zj  is  the  measurement  at  the  reference  sphere,  Cat  is  the  correction  factor  for  the 
sample  sphere  measured  with  no  sample  at  port  I,  port  II  covered  and  the  laser  off,  while 
Czr  is  the  correction  for  the  light  in  the  reference  sphere  with  the  laser  off.  The  2.33  term 
is  70/30  to  account  for  the  beam  splitter  difference. 


Pig-tailed  Diode  Laser 


Oscilloscope 


Fig.  1.  Setup  for  diffused  reflectance  and  total  transmittance  measurements. 
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Glass  Window 


Heater  Elements 


Fig.  2.  Temperature  controlled  sample  holder  for  diffused  reflectance  and  total 


Inverse  Adding  Doubling  (IAD)  Method 

The  IAD  computer  program  was  used  to  solve  the  radiative  transport  equation. 
The  IAD  computer  program  must  be  supplied  with  the  experimental  values  of  total 
diffuse  reflectance  (f?d)  and  total  transmittance  (Tt)  along  with  values  of  index  of 
refraction  (n)  and  scattering  anisotropy  (g)  of  the  sample.  The  IAD  program  guesses  a  set 
of  optical  properties  ( a  and  t,  defined  below)  and  then  calculates  values  for  /?,i  and  Tt. 
This  process  is  repeated  until  the  calculated  and  measured  values  of  Rd  and  Tt  are  within  a 
specified  tolerance.  Values  of  a  and  r  provided  by  the  IAD  method  are  used  to  calculate 
the  absorption  coefficient  (//a)  and  scattering  coefficient  (//s)  for  individual  tissue  using 
Eqs.  (3)  and  (4). 


a  =  Us/ (  jUs  +  JU a) 

(3) 

T=t(jUs  +  jU  a), 

(4) 

In  the  above  t  is  the  physical  thickness  of  the  sample  measured  in  cm. 

According  to  Prahl  et  al.1,  the  validity  of  the  IAD  method  for  samples  with 
comparable  absorption  and  scattering  coefficients  is  especially  important,  since  other 
methods  based  on  only  diffusion  approximation  are  inadequate.  They  have  also  stated 
that  since  both  anisotropy  phase  function  and  Fresnel  reflection  at  boundaries  is 
accurately  approximated,  the  IAD  method  is  well  suited  to  optical  measurements  for 
biological  materials  sandwiched  between  two  glass  slides.  Details  of  the  IAD  method 
have  been  previously  provided  by  Prahl  et  al. 
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Sample  Preparation 

One  pair  of  bovine  eyes  was  procured  from  a  local  slaughterhouse.  This  pair  of 
eyes  was  enucleated  fresh,  packed  on  ice  and  immediately  transported  to  the  research 
laboratory  (within  one  hour).  Immediately  upon  arrival  at  the  research  laboratory,  fresh 
(not  frozen)  whole  eyes  were  dissected  on  ice  to  remove  vitreous  humor  for  optical 
measurements.  Approximately  three  milliliters  of  vitreous  humor  was  removed  from  the 
bisected  bovine  eye  using  a  Ice  tuberculin  syringe,  sans  the  fine  needle.  Vitreous  humor 
was  then  introduced  into  the  sample  holder  described  above  and  shown  in  Fig.  2.  Visual 
inspection  showed  clear  vitreous  humor  at  the  time  when  optical  measurements  were 
conducted. 

Experimental  Results 

The  diffuse  reflectance  (Rd)  and  total  transmittance  (7i)  of  the  bovine  vitreous 
humor  at  980  in  the  temperature  range  of  24  to  40°C  are  given  in  Figs.  3  and  4, 
respectively.  There  seemed  to  be  little  correlation  between  Rd  or  Tt  and  temperature.  A 
repeat  of  the  total  transmittance  measurements  at  23°C  at  five-minute  intervals,  roughly 
the  same  duration  between  temperature  readings,  for  ten  readings  is  shown  in  Fig.  4. 
Note  the  spread  of  the  data  at  this  one  temperature  is  close  to  that  of  the  data  taken  at 
other  temperatures.  It  was  thus  concluded  that  any  temperature  dependence  of  the  diffuse 
reflectance  and  total  transmittance  could  not  be  discerned  with  the  current  setup.  No 
attempt  was  made  at  any  other  wavelengths  for  any  sample. 

The  absorption  and  scattering  coefficient  for  the  bovine  vitreous  humor  from  the 
data  averaged  over  all  temperatures  was  determined  to  be  0.75  cm'1  and  0.03  cm"1, 
respectively.  The  refractive  index  and  anisotropy  were  considered  to  be  independent  of 
temperature  and  taken  from  Sardar  et  al."  to  be  1.33  and  0.99,  respectively.  Values  for 
absorption  and  scattering  coefficient  for  bovine  vitreous  humor  at  room  temperature 
reported  by  Sardar  et  al.  are  1.27  cm"1  and  0.88  cm'1,  respectively.  Though  the  absorption 
coefficients  are  comparable  the  scattering  coefficients  are  not.  It  is  possible  that  some 
specular  reflection  was  inadvertently  measured  in  the  work  done  by  Sardar  et  al. 
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Temperature  (  C) 


Fig.  4.  Total  transmittance  (Tt)  data  for  bovine  vitreous  humor  at  980nm.  The  dots  are 
the  average  and  the  bars  are  the  standard  deviation  of  10,000  data  points. 


Temperature  Dependence  of  — — 

dT 


Experimental  Details 

In  the  current  work,  we  use  a  Mach-Zehnder  interferometer  to  measure  the 

c)ti 

—  (T )  parameter  without  knowing  the  actual  value  of  the  refractive  index.  The 
dT 

schematic  of  which  is  shown  in  Fig.  5  where  BS1  and  BS2  are  beam  splitters  and  Ml  and 
M2  are  mirrors.  The  paths  BS1-M1-BS2  and  BS1-M2-BS2  are  referred  to  as  the 
reference  arm  and  the  sample  arm  of  the  interferometer,  respectively. 

The  tissue  was  placed  in  the  sample  arm  of  the  interferometer  using  a 
temperature-controlled  sample  holder,  the  details  of  which  are  shown  in  Fig.  6.  The 
sample  holder  was  constructed  from  aluminum  9.67  mm  thick  with  an  aperture  of  about 
12  mm  in  diameter  and  had  two  1  mm  thick  glass  windows.  Two  electric  heating 
elements  (Omega  SRFG101)  were  affixed  to  either  side  of  the  sample  holder.  The  current 
through  the  heater,  and  thus  the  heat  introduced  into  the  system,  was  controlled  by  a 
voltage  regulator  (Powerstat  3PN1 16B).  The  temperature  of  the  sample  was  monitored  by 
an  electronic  thermometer  system  (Flart  Scientific  1006).  In  the  current  study  data  was 
recorded  from  37  to  41°C  at  0.05°C  intervals.  The  interferograms  were  captured  by  a 
CCD  camera  (Spiricon  Cohu  4800)  and  recorded  as  bitmap  images  for  final  analysis.  A 
typical  fringe  interferogram  for  the  setup  used  here  is  shown  Fig.  7.  The  light  source  of 
the  interferometer  was  a  5  mW  helium-neon  laser  (Electro-Optics  LFIRP-0501)  at 
wavelength  632.8  nm. 


Laser 


BS1 


Macro  Lens 
CCD  Camera 


Fig.  5.  Mach-Zehnder  interferometer  setup  for  measuring 
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Heater  Elements 


Fig.  6.  Temperature  controlled  sample  holder  for  — ( T )  measurements. 

dT 


Fig.  7.  Typical  fringes  observed  with  the  setup  in  Fig.  5. 
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Principle 

For  our  experimental  conditions,  the  thickness  of  the  sample  in  the  optical 
pathway  was  9.67  mm,  and  we  assumed  that  the  sample  difference  in  the  direction  of  the 
optical  pathway  is  negligible.  In  such  case  we  can  treat  the  problem  two  dimensionally  so 
that  every  point  in  the  observed  region  can  be  expressed  by  an  orthogonal  coordinate 
(x,y).  If  we  study  the  sample  only,  and  suppose  that  the  initial  refractive  index  of  the 
solution  is  n(x,y,T0),  after  a  small  change  in  temperature  AT  the  refractive  index 

changes  to  n(x,y,T0  +  AT )  at  the  point  (x,y).  The  phase  change  at  point  (x,y)  is  thus 
<j)(x,yj0+  AT)- (j)(x,  y,T0),  and  we  have: 


[n(x,  y ,  T0  +  AT)  -  n(x,  y ,  T0  )]d  =  [(/){x,  y,  T0+AT)-  (f>(x,  y,  T0 )]  (5) 

2n 

Dividing  both  sides  of  equation  5  by  the  thickness,  d,  and  AT  and  taking  the  limit  as  AT 
tends  to  0,  we  get: 


dn 

dT 


(■ x,y,T ) 


X  8<j> 
2nd  dT 


{x,  y,  T ) 


(6) 


Fixing  the  point  (x,y) 


we  get: 
dn 


dT 


X  d<t>xy 
2nd  dT 


(r) 


(7) 


Hence,  once  —  (r)  is  known,  —  (t)  can  be  determined. 
dT w  dTK 

—  (T )  from  an  interferometric  setup  as  discussed  below. 
dT 


We  determine  the  tenn 


For  a  50-50  interferometer,  the  interference  tenn  for  the  intensity,  Tnt,  is  given  by: 

7„u  =  21  cos(d</>)  (8) 

where  /  is  the  laser  intensity  and  5(j)  is  the  phase  difference  between  the  two  paths.  For 
wide  field,  temperature  dependent,  equation  8  takes  the  fonn: 

4t  4,  y,  T)  =  21  cos  (s<t(x,  y,  T ))  (9) 

Fixing  the  point  (x,y)  we  can  expand  the  phase  difference  as: 

S</>x  V(T)  =  a' o+aj1  +  a2T2  +  a3T3  (10) 


8 


Since  the  reference  arm  of  the  interferometer  is  at  a  fixed  temperature,  we  can  absorb  the 
difference  in  phase  at  the  initial  temperature  into  the  constant  a\  to  get  another  constant 

a0 ,  and  thus  we  can  recast  equation  10  as: 

<f>x  v  (71)  =  <70  +  T  +  a -,T ~  +  g2T2  (1 1) 


where  (j)x  v  ( T )  is  now  the  phase  of  the  light  passing  through  the  sample  relative  to  the 
reference  path  as  a  function  of  sample  temperature. 

Thus  for  a  given  point  (x,y)  we  have  the  model: 

Iiatxy(T)  =  max(o,  A  cos(a0  +  axT  +  a2T2  +  a3T3)+  b)  (12) 


where  the  B  term  is  needed  to  account  for  the  threshold  to  be  overcome  in  the  detector 
and  the  max  term  is  needed  since  the  intensity  cannot  be  less  than  zero.  It  is  this  model 
that  we  fit  observed  data  from  a  Mach-Zehnder  interferometer  to  obtain  <f>(T)  which  is 


then  differentiated  to  obtain  — M7T  then  scaled  by - 

8T  2nd 


(equation  7)  to  obtain  ~{T )  • 


Sample  Preparation 

Two  pairs  of  bovine  eyes  were  procured  from  a  local  slaughterhouse.  These  pairs 
of  eyes  were  enucleated  fresh,  packed  on  ice  and  immediately  transported  to  the  research 
laboratory  (within  one  hour).  Immediately  upon  arrival  at  the  research  laboratory,  fresh 
(not  frozen)  whole  eyes  were  dissected  on  ice  to  remove  aqueous  vitreous  humor  for 
optical  measurements.  Approximately  three  milliliters  of  aqueous  humor  was  removed 
from  the  two  bisected  bovine  eye  using  a  lcc  tuberculin  syringe.  Aqueous  humor  was 
then  introduced  into  the  sample  holder  described  above  and  shown  in  Fig.  6.  For  the 
vitreous  humor,  approximately  three  milliliters  of  vitreous  humor  was  removed  from  a 
bisected  bovine  eye  using  a  lcc  tuberculin  syringe,  sans  the  fine  needle.  Vitreous  humor 
was  then  introduced  into  the  sample  holder  after  the  experiment  was  perfonned  for  the 
aqueous  humor  and  the  sample  holder  washed.  Visual  inspection  showed  clear  aqueous 
and  vitreous  humor  at  the  time  when  optical  measurements  were  conducted. 

System  Calibration 

As  a  means  of  calibrating  our  setup,  we  relied  on  the  data  of  Schiebener  et  al.3  To 
obtain  a  closed  form  solution  of  the  refractive  index  of  water  at  632.8  nm  at  atmospheric 
pressure  we  best  fit  the  data  reported  by  Schiebener  et  al.  to  polynomials  of  second  and 
third  order.  Both  polynomials  give  a  good  fit  to  the  refractive  index  data  for  water 
reported  in  Schiebener  et  al.  as  shown  in  Fig.  8.  Differentiating  the  best- fit  polynomials 
however  reveals  considerable  difference  between  the  second  and  third  order  polynomials 
as  shown  in  Fig.  9.  Significant  difference  between  the  two  occurs  in  our  temperature 
range  of  interest  of  37  to  41°C  (see  Fig.  10).  We  hence  use  a  third  order  polynomial  as 
the  true  fit  for  the  data  reported  in  Schiebener  et  al. 
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Fig.  8.  Refractive  index  values  for  water  from  Schiebener  et  al.  and  best  fit  second  and 
third  order  polynomials. 


Fig.  9.  Computed  — ( T )  curves  for  the  best  fit  second  and  third  order 
dT 

polynomials  to  Schiebener  et  al.3  from  0  to  100°C. 
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Fig.  10.  Computed  — ( T )  curves  for  the  best  fit  second  and  third  order 
dT 

polynomials  to  Schiebener  et  al.3  from  37  to  41°C. 


Experimental  Results 

For  all  experiments  the  voltage  across  the  four  resistive  heaters  was  set  to  30.0 
volts  resulting  in  a  temperature  increase  rate  of  about  0.005°C/second  throughout  the 
temperature  range  of  37  to  41°C.  The  intensity  captured  by  the  CCD  camera  for  water 
for  a  given  pixel  is  shown  in  Fig.  1 1  as  solid  dots.  The  best  fit  to  this  curve,  obtained  by 
a  genetic  algorithm  written  in  Matlab,  for  equation  13  is  shown  as  the  solid  line.  The 
difference  in  the  calculated  and  the  fitted  data  from  Schiebener  et  al.3  is  shown  as  the 
green  line  is  Fig.  12.  This  is  taken  as  the  correction  factor  for  our  experimental  setup, 
which  is  thought  to  be  due  to  the  phase  shift  induced  by  the  glass  windows  on  either  side 

c)ti 

of  the  sample.  A  comparison  to  literature4  gives  — ( T )  for  glass  to  be  of  the  order  of 

dT 

magnitude  of  10'5,  comparable  to  our  results.  The  same  experiment  for  bovine  aqueous 
humor  and  vitreous  humor  at  632.8  nm  were  performed.  Figs.  13  and  15  show  the 
experimental  values  and  best  fit  for  the  intensity  as  a  function  of  temperature  for  bovine 
aqueous  humor  and  vitreous  humor,  respectively.  Figs.  14  and  16  show  the  calculated 

dn 

—  ( T )  for  the  measured  samples  plus  the  corrected  values  for  bovine  aqueous  humor 
dT 

and  vitreous  humor,  respectively. 
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For  bovine  aqueous  humor  the  best  fit  to  measured  data  yielded: 


——(T )  =  8.268 x  10  4  -  4.313  x  +  4.960 x  10  7 T2  (13) 

dT 


For  bovine  vitreous  humor  the  best  fit  to  measured  data  yielded: 

—  (r)  =  5. 815x10  3  -  3.000  xl0~4r  + 3.758x10 ~6T2  (14) 

dT 


These  values  are  of  the  same  order  of  magnitude  of  water  as  shown  in  Fig.  17.  Note  that 

c)ti 

—  (T )  for  bovine  aqueous  humor  is  close  to  that  of  water.  Depending  on  the  precision 
dT 

required  for  the  value  of  — (T )  it  may  suffice  to  take  — (T )  for  bovine  aqueous  humor 

dT  dT 

to  be  that  of  water. 


Fig.  11.  Data  and  best  fit  for  interference  of  water. 
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Fig.  13.  Data  and  fit  for  bovine  aqueous  humor  at  632.8  mn. 
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Fig.  14.  Measured  and  corrected  — ( T )  for  bovine  aqueous  humor  at  632.8  nm. 

dT 
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Fig.  15.  Data  and  fit  for  bovine  vitreous  humor  at  632.8  nm. 
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Fig.  16.  Measured  and  corrected  — ( T )  for  bovine  vitreous  humor  at  632.8  nm. 

dT 


Fig.  17.  Measured — -( T )  for  water  and  bovine  aqueous  and  vitreous  humor  at  632.8  nm. 
dT 
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Conclusion 

The  objective  of  this  research  was  to  systematically  measure  the  temperature 
dependence  of  optical  properties  of  various  bovine  eye  tissues  using  different 
experimental  and  computational  techniques  at  elevated  temperatures  (above  body 
temperature). 

The  diffuse  reflectance  (Rd)  and  total  transmittance  (7))  of  the  bovine  vitreous 
humor  at  980  in  the  temperature  range  of  24  to  40°C  were  measured.  It  was  determined 
that  any  temperature  dependence  of  the  diffuse  reflectance  and  total  transmittance  could 
not  be  discerned  with  the  current  setup  in  the  temperature  range  of  interest.  The  average 
Rd  and  T{  values  over  all  measured  temperatures  were  input  into  the  inverse  adding 
doubling  (IAD)  computer  program  yielding  absorption  and  scattering  coefficients  of  0.75 
cm’1  and  0.03  cm’1,  respectively.  Values  for  absorption  and  scattering  coefficient  for 
bovine  vitreous  humor  at  room  temperature  reported  by  Sardar  et  al.  are  1 .27  cm'  and 
0.88  cm’1,  respectively.  Though  the  absorption  coefficients  are  comparable  the  scattering 
coefficients  are  not.  It  is  possible  that  some  specular  reflection  was  inadvertently 
measured  in  the  work  done  by  Sardar  et  al.  No  attempt  was  made  to  measure  temperature 
dependent  absorption  and  scattering  coefficients  at  any  other  wavelengths  for  any  sample. 

QH 

The  rate  of  change  of  the  refractive  index  with  respect  to  temperature,  — ,  was 

dT 

performed  at  632.8  nm  for  bovine  aqueous  and  vitreous  humor  in  the  temperature  range 
of  37  to  41°C  using  a  Mach-Zehnder  interferometer  setup.  The  system  was  calibrated  for 

water  using  data  published  by  Schiebener  et  al.3  The  measured  — (7)  for  bovine 

dT 

aqueous  humor  and  vitreous  humor  were  of  the  same  order  of  magnitude  of  water.  The 

dn 

measured  — ( T )  of  bovine  aqueous  humor  was  comparable  to  that  of  water,  while  the 
dT 

dn 

measured  — ( T )  for  bovine  vitreous  humor  was  slightly  less.  These  values,  as  vet  not 
dT 

reported  in  literature,  are  necessaty  for  accurately  modeling  thermal  lensing  and  hence 
light  propagation  in  the  eye. 
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